Singular Scaling Functions in Clustering Phenomena 
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We study clustering in a stochastic system of particles sliding down a fluctuating surface in one 
and two dimensions. In steady state, the density-density correlation function is a scaling function of 
separation and system size. This scaling function is singular for small argument — it exhibits a cusp 
singularity for particles with mutual exclusion, and a divergence for noninteracting particles. The 
steady state is characterized by giant fluctuations which do not damp down in the thermodynamic 
limit. The autocorrelation function is a singular scaling function of time and system size. The 
scaling properties are surprisingly similar to those for particles moving in a quenched disordered 
environment that results if the surface is frozen. 
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I. INTRODUCTION 

Clustering phenomena are widespread. On the one 
hand, they include the growth of regions of a homoge- 
neous ordered phase, such as droplets of a hquid or do- 
mains in a magnet; on the other hand, they also include 
the formation of inhomogeneous structures such as those 
that arise in granular gases. In general, clusters break 
and re-form due to the action of noise, whether thermal 
or induced by external driving. Concomitantly, the dis- 
tribution of clusters - their sizes and spatial locations - 
keep evolving in time. A useful diagnostic of the nature of 
clustering, both in steady state and during the approach 
to it, is provided by the two-point density-density corre- 
lation function. In this paper, we will discuss the scaling 
properties of correlation functions, and show that clus- 
tered states of different sorts correspond to different types 
of singularities in the corresponding scaling functions. 

One of the most familiar and well studied examples of 
clustering is afforded by phase ordering kinetics. Here 
the clusters are simply domains of ordered phases which 
form following a rapid quench from an initially disordered 
state to a regime which favours ordering. The coarsening 
of ordered domains in time is captured by the two-point 
correlation function; it shows simple scaling properties 
when spatial separations are scaled bythe mean domain 
size, itself a growing function of time [l[ . In steady state, 
the system size replaces the domain size in the scaling 
function. 

There are a number of situations in which particles 
cluster in such a way that the resulting density profile is 
patchy and far from homogeneous. For instance, parti- 
cles may arrange themselves in a critical state with power 
law decays of the correlation function, e.g. molecules in a 
fluid at the critical point [3] or galaxies in the universe Q . 
However, the sort of clustering of interest in this paper 
is different in that it is associated with long range order, 
unlike a critical state, yet has strong fluctuating prop- 
erties which make it quite different from a conventional 
phase ordering system. 

We focus on particles subjected to randomly fluctu- 



ating forces, which result in an inhomogeneously clus- 
tered state. Particles advected by fluid flows exhibit a 
clustering tendency whose strength depends on the com- 
pressibility of the fluid and the inertia of the particles 
til tSi, ^] • Earlier studies have shown that in the absence 
of noise, particle trajectories can coalesce as time ad- 
vances [3, Q . Not surprisingly, external noise — which 
allows coalesced clusters to break up — has a profound 
effect, and the dynamical steady state which results from 
the balance between making and breaking clusters has 
interesting characteristics. Our aim in this paper is to 
discuss the nature of this state, with a focus on the scal- 
ing properties of correlation functions. We study simple 
stochastic models of particles sliding down a fluctuating 
surface H, [M [H [H [H, [II • As for the case of phase 
ordering kinetics, the two-point correlation function has 
simple scaling properties. However, now the associated 
scaling function develops singularities - these are cusp 
singularities for particles with mutual exclusion, and di- 
vergences for noninteracting particles. The occurrence of 
such singularities is quite robust, and survives as param- 
eters of the model and the dimension are changed. More- 
over, recent work has shown that such singularities are 
found also in other completely different systems which 
exhibit inhomogeneous clustering [l^, [l^ . 



In a finite system, a dynamical steady state is reached 
in which correlations are scaling functions of the separa- 
tion and system size. Like conventional phase ordering 
systems, this system exhibits long range order, but un- 
like them, the scaling functions are singular at small ar- 
gument, indicative of a steady state which is intrinsically 
different from a phase separated state in an equilibrium 
system. A crucial point about the steady state is the 
existence of strong fluctuations which do not die down 
in the thermodynamic limit, due to which the degree of 
ordering in the system keeps fluctuating, but never van- 
ishes. 
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II. SCALING IN PHASE ORDERING SYSTEMS 

Phase ordering kinetics deals with the onset of order in 
a system which is brought rapidly from a disordered state 
to an order-promoting regime [l| . Domains of competing 
phases coarsen in time, a process which is captured by the 
two-point correlation function C{r,t) = {(Ji(t)ai+r{'t))) 
where (7^ = ±1 is an Ising variable which describes a 
spin in a magnet, or an occupancy variable (1 + ai)/2 
within a lattice gas description. In an infinite system, 
the correlation function has the scaling form 

Csir,t)^Y{r/C{t)) (1) 

in the scaling limit r — > oo, t — > cxd with the ratio y = 
r/C{t) held fixed. Here C{t) is a growing, time-dependent 
length scale which describes the typical linear size of a 
cluster at time t] typically it grows as a power law in 
time, C{t) - tl^' . 

The scaling form of Eq. (1) holds only if the separa- 
tion r 3> the correlation length ^. In general, the cor- 
relation function C(r, t) contains both the scaling part 
Cs{r,t) and an analytic part (^^(r) which comes from 
short-distance (r ^ ^) correlations, independent oi C{t). 
In practice, it is important to account for the occurrence 
of Ca {t) while analyzing data for the correlation function 
to see whether scaling holds. 

When the system size L is infinite, the domain size 
C{t) can grow indefinitely, with the scaling form of Eq. 
1 continuing to hold. However, if L is finite, then Eq. 
1 ceases to hold once t is large enough that C{t) « L. 
Beyond this, the system reaches steady state and Eq. 
(1) is replaced by 

Cs{r,L) = Y{r/L). (2) 

There are two important characteristics of the scaling 
functions Y{y) of Eq. 1 and Y{y) of Eq. 2: (i) the scal- 
ing function has a finite intercept, (ii) the scaling function 
falls linearly for small y. Property (i) is important as the 
intercept is a measure of long range order (LRO). To see 
this for Y{y), recall that LRO is defined in an infinite 
system in steady state as mp = Yvctir ^ oo {cf iU i+r) ■ Now 
consider two points separated by an arbitrarily large but 
fixed distance r in a coarsening system. As t — > cx), do- 
main sizes grow without bound and the interior of each 
domain is well approximated by steady state. Thus at 
large times, both points would belong to the same do- 
main with probability one, implying that their correla- 
tion equals to§. Since t ^ oo leads to ?/ ^ 0, no matter 
how large r is, this implies that the intercept in Fig. 1 is 
ml. 

Further, property (ii) which describes the manner in 
which the scaling function falls from the intercept value 
mQ is also significant. For scalar order parameters, it 
falls linearly: Y{y) « ml — b\y\ as y For 
a d-dimensional system, this translates into a decay ~ 
[fc£(i)]~^'^+^' for the structure factor at large wave- vector 
k. This power law decay is known as the Porod Law p7j . 




(a) (b) (c) 

FIG. 1: The behaviour of the scaling part of the correlation 
function is depicted for (a) normal phase ordering systems (b) 
a system of particles with exclusion, sliding down a fluctuating 
surface (Model E) (c) noninteracting particles sliding down a 
fluctuating surface (Model F). The linear drop of Cs in (a) 
[the Porod law] is replaced by a cusp singularity in (b), and 
by a divergence in (c). 



The linear fall of correlations for \y\ = \r/C\ ^ 1 can be 
traced to the existence of well-defined interfaces between 
two phases: Two points a distance r apart are positively 
correlated if both belong to the same domain, and neg- 
atively correlated if a single interface intersects the line 
joining them; for small values of r/£, the chance of this 
happening is proportional to r and inversely proportional 
to C, accounting for the linear drop of the scaling func- 
tion at small argument. 

In subsequent sections we investigate how far these 
properties are respected when we have clustering of an 
inhomogeneous sort. We find that there is a class of inho- 
mogeneous clustering phenomena, exemplified by models 
of particles sliding on randomly fluctuating surfaces, for 
which most of the properties discussed above are valid, 
except for the form of singular behaviour of the scaling 
function at small argument. These models and a discus- 
sion of the singularities they display are discussed in the 
following sections. 



III. MODEL: PARTICLES SLIDING ON A 
FLUCTUATING SURFACE 

The model consists of a system of overdamped parti- 
cles sliding down under gravity along the local slope of a 
fluctuating interface (Fig. 2) [all^- In 1-d, we model the 
interface through the single-step model J^. It consists 
of a series of links, with the slope of the link between the 
i'th and {i + l)'th site being T^^l = ±1, so that a typical 
{t} configuration is ///\\\\\/\/- The height at site i 
is given by hi — The elementary dynami- 

cal move is a stochastic corner flip, which involves the 
exchange of adjacent r's: the transition /\ \/ occurs 
with rate pi, while the reverse transition \/ /\ occurs 
with rate qi. For large r and t, with pi = qi = 1 the inter- 
face dynamics reduces to the Edwards- Wilkinson model 
, while with pi ^ oj, it reduces to the Kardar-Parisi- 
Zhang (KPZ) model [20|. Note that the overall slope of 
the interface T = X^i^i ''"i-i-J" conserved; we use peri- 
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odic boundary conditions and work with interfaces with 
no overaU tilt i.e. T — Q. 

For the dynamics of the particles, there are two cases: 
particles with mutual exclusi on ( E) [9| , and particles with 
no interactions, i.e. free (F) p^ . 



(a) \/ 



O . O /V >o 

(b) \:/ 



FIG. 2; Typical configurations and moves for lattice mod- 
els of particles sliding down a fluctuating surface in l-d. (a) 
corresponds to Model E (particles with exclusion) and (b) to 
Model F (particles with no interactions). 

In case E (Fig. 2a), a particle and hole at sites i and 
i + 1 exchange at a rate which is governed by the slope 
variable Tj_|_i. A move down the slope (•\o o\» or 
o/» — > •/o) occurs with rate p2, while the reverse transi- 
tions occur with rate q2 < P2- 

In case F (Fig. 2b), there is no limit on the occupancy 
of any site, and a particle moves to its neighbouring site 
with rate p2 or q2 depending on whether it traverses the 
intervening t bond in the downward or upward direction. 
If q2 = 0, particles move only downward and coalesce at 
the bottoms of valleys. Nevertheless, particle clusters can 
break again if they are carried, by surface evolution, to a 
local hilltop; subsequently, particles are equally likely to 
fall left or right, leading to a break-up of the cluster. 

Notice that in both cases E and F, the particles are pas- 
sive scalars: while the evolution of particles occupancies 
{rii} is governed by the interface configuration {t^^i}, 
the evolution of the interface is autonomous and is not 
influenced by the particles. 

Model E with EW dynamics is related to the Lahiri- 
Ramaswamy (LR) model of sedimenting colloidal crystals 
(21I |2^ . In the present context, the general LR model 
allows for the particles to act back on the fluctuating 
surface, in which case phase-separated and disordered 
phases ensue [12, HI) . The case of passive particles under 
consideration lies at the boundary between these more 
conventional phases [23 |. 

A lattice model closely related to model F with KPZ 
dynamics was studied by Drossel and Kadar [2^, \2^ : 
their results are discussed in Section 5 along with ours. 
Interestingly, model F with KPZ dynamics relates to the 
problem of passive scalar advection by a noisy Burgers 



flow. This follows from the well known mapping [20] be- 
tween the KPZ equation and the noisy Burgers equation 
which describes a compressible flow, and from the obser- 
vation that movement of a particle along the local slope 
in the present context translates, through the mapping, 
to a particle being advected by the local velocity field. 



IV. SLIDING PARTICLES WITH EXCLUSION 

We turn to the properties of a system of particles with 
mutual exclusion on a fluctuating interface, evolving un- 
der the dynamics discussed in the previous section. 



A. Coarsening 

At t = 0, we start out with a configuration of the inter- 
face drawn from steady state, and particles distributed 
randomly on sites. Numerical simulations of this process 
using a very large system size were performed to monitor 
the two-point correlation functions C (r, t) , averaged over 
many histories [1, [l^]. Particles and surface elements 
were updated at equal rates. The results show that there 
is a good scaling collapse of the data to the form of Eq. 
1, with jC{t) ^ t^/^ where z is the dynamical exponent 
which characterizes the interface dynamics. This is be- 
cause in time t, new valleys of size t^/^ form, so that 
particles in basins of that size are drawn together. In 1 
dimension, we have z — 2 for EW evolution, and z — 3/2 
for KPZ evolution. 

The significant new point is that the scaling function 
has a cusp singularity 



Cs{y) ^ ml ^ b\y\" asy^O 



(3) 



with y ~ r/C{t). This translates into a tail ^ {kC{t)Y^°' 
for the wave- vector dependence of the structure factor, 
implying non-Porod behaviour. The cusp exponent a is 
found to be ~ 0.5 with EW evolution and c± 0.25 with 
KPZ evolution of the interface The value of a is 

found to be unaffected by variations of the density of 
particles or the relative ratio of interface and particle 
update rates. Moreover, the phenomenon persists with 
2-d interfaces as well [llj; with KPZ dynamics, it is found 
that a ~ 0.38 (see Section 4.4). 



B. Steady state and Cluster size distribution 

A study of finite systems of size L shows that in steady 
state, the correlation function is a scaling function of sep- 
aration and size, as in Eq. 2 fTo'|. Once again, the scaling 
function shows a cusp singularity, and the value of the 
cusp exponent a matches that obtained in the coarsening 
studies discussed in Section 4.1. 

A clue to the nature of the state comes from a study 
of the distribution P{i) of cluster sizes I A clus- 

ter is defined as a stretch of continuously occupied sites. 
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with unoccupied perimeter sites. In the 1-d EW case, it 
is found that P{i) ~ with 9 ~ 1.8; symmetry im- 
plies the same form for hole clusters. In the 1-d KPZ 
case, simulations reveal that hole clusters follow a power 
law decay with 6 ~ 1.85, whereas the decay for particle 
clusters do not seem to follow a simple power law. 

The power-law decay of P{£) can be related to the cusp 
observed in the two-point correlation function |1Q] within 
the independent interval approximation (llA) which is 
based on the premise that the occurrence of successive 
intervals (clusters of particles and holes) are independent 
events. With this approximation, the Laplace transforms 
Pis) = J^M e-^'P{t) and C{s) = dr C{r) are re- 
lated by p3| 



1 



sC{s 



2 1 - P(s) 
<l>l + P{s) 



(4) 



where < £ > is the mean cluster size. For the slow power 
law decay P{t) ~ t^^ , where £ is limited by £ = i, we 
have < ^ >« aL2-» and P(s) w 1 - hs^-^ for l/L < 
s <C 1. Equation 4 then leads to the scaling from of Eqs. 
2 and 3, with the cusp exponent given by 



a = 2 - I 



(5) 



The occurrence of LRO in this fluctuating system can 
also be understood by noting that extremal statistics ap- 
plied to P{t) ^ with \ < 6 <2 leads to the conclu- 
sion that the largest clusters are of size ~ L. 



C. Adiabatic limit: Coarse-grained depth (CD) 
model 

We can calculate 9 and thus a within the IIA for a 
class of models motivated by considering the adiabatic 
limit of infinitely slow interface evolution. Then the 
height field {hi} is a quenched random variable with 
hi = '^i-j- '^■'^^ steady state of the particle distri- 

bution is a thermal equilibrium state with temperature 
defined by T = V£n{p2/q2) where V is the potential en- 
ergy drop across a single link. In the limit of low T, this 
equilibrium state is characterized by a filling level akin 
to the Fermi level, namely the height up to which the 
particles are filled. In this limit, the correlation function 
< SjSj+r > coincides with C(r). 

This motivates us to define and study more general 
coarse-gained depth models (CD models) [S] with vari- 
able filling levels 



i{t) = sgn [hi{t) - h] 



(6) 



where h is the height of a surface cut. Another natural 
choice for h is the instantaneous average height of the 
configuration {hi}, but with this definition the cut need 
not coincide with the filling level defined in the previous 
paragraph. Yet another possibility for the cut is to pin it 
to the instantaneous height of the first site. For the 1-d 



EW and KPZ interfaces, the configuration of the inter- 
face is given by the trajectory of a random walk and the 
filling level defines a cut of this trajectory. The intervals 
between successive crossings then follow a distribution 
p{i) ^ £~^^'^. Using Eq. 5, this leads to a cusp exponent 
a ^ 1/2. 



D. Higher dimensions 

Singularities in scaling functions are also found to oc- 
cur in a system of particles slidin g d own a fluctuating 
rough surface in two dimensions [ll|. The surface is 
modeled by a discrete solid-on-solid model with a sim- 
ple growth rule: a site on a square lattice is selected a 
random and its height is decreased by 2 units, provided 
the height at all four of its neighbouring sites is lower. 
The asymptotic properties of this model are expected 
to be the same as those of the (2-|-l)-dimensional KPZ 
equation. Particles are initially distributed at random, 
and updated by choosing a random particle, attempting 
to move it to a randomly chosen neighbour, and actually 
moving it provided the local slope is downward and the 
target site is unoccupied. 

Results on coarsening [llj indicate that Eq. 1 holds, 
with £(t) ~ i^/^ and z ~ 1.6, close to the value of the dy- 
namical exponent which defines the relaxation of a (2-1-1)- 
dimensional KPZ surface. The cusp exponent a char- 
acterizing the decay of the scaling form of the density- 
density correlation function is ~ 0.38. A study of the 
CD model (with h taken to be the average height) shows 
that Eq. 2 holds, with a ~ 0.43. 

In steady state, we may relate the cusp exponent a in 
the CD model to first return probabilities through the fol- 
lowing argument [Tl| . Consider a self-affine surface with 
roughness exponent x, and define P{£) as the probability 
that the surface first returns to its starting height a dis- 
tance £ along an arbitrary linear direction. P{£) decays 
as power law for £ < L : P{£) ~ ^"(^-x). Qn applying 
the IIA to different segments on a straight line connect- 
ing two points in a distance r apart, Eq. 5 leads to the 
conclusion that a = X: i-e. the cusp exponent for the 
CD model is determined by the roughness exponent of 
the fluctuating, driving surface. This is borne out by the 
results in both 1 and 2 dimensions (a = x = 1/2 in 1-d; 
a ~ 0.43, X =i 0-4 in 2-d). 



E. Dynamics 

As with static (equal-time) correlation functions in 
steady state, dynamical properties also exhibit scaling 
characterized by singular scaling functions 14]. The 
autocorrelation function A(t) — {<Ji(0)ai{t)) was moni- 
tored by Monte Carlo simulation, and found to be de- 
scribed by a scaling function W{t/L^), in the scaling 
limit t oo,L — > oo with i/L^ held finite. Parallel to 
the static situation, the scaling function Ws{w) exhibits 
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a cusp singularity at small argument: 

Ws{w) Kml{l-h'\wf) 



w 







(7) 



where w = t/ , ttiq measures long-range order and (3 is 
the cusp exponent. The results of simulations indicate 
that (3 ~ 0.22 for Edwards- Wilkinson surface dynamics 
and P ~ 0.18 for KPZ surface dynamics [T^ . 

Just as for static correlations (section 4.3), insight 
into singular scaling for dynamic correlations can be ob- 
tained from analytic calculations on the corresponding 
CD model. To this end, the autocorrelation function 
(si(O)si(t]) of the CD model with EW evolution was cal- 
culated [iJ] and shown to follow the form of Eq. 2 with 
mo = 1 and /3 = 0.25. 

Further, aging functions of the form A{ti,t2) = 
{(Ji{ti)ui{t2)) were investigated for a coarsening system, 
starting from a state with randomly distributed parti- 
cles [lj|. If ti,t2 ^ 1, the aging function ^ is a func- 
tion of the ratio ti/t2 alone. For ti t2, we find 
A ~ TOg[l — 6i(i2/ii)'^], whereas for ti <^t2, the function 
A follows a power law decay {ti/t2)^ ■ Simulations show 
that 7 ~ 0.69 for EW interface dynamics, 7 ~ 0.82 for 
KPZ dynamics, and 7 ~ 0.82 for the CD model. 



F. Ordering and Giant fluctuations 

As we have seen, the two-point correlation function 
characterizing the steady state approached by a system 
of particles driven by a fluctuating surface exhibits long 
range order. Let us ask for a one-point function or order 
parameter which can be used to characterize the steady 
state. To this end, we monitored the Fourier transform 



Q{k) = 



1 ^ 



2nm 



(8) 



where m — 1, 2, • • • , L — 1. The expectation value Q* — 
(Qi) of the longest-wavelength mode Qi = Q{2n/ L) cor- 
responding to TO = 1 is a putative order parameter. For a 
state which exhibits complete phase separation in a half- 
filled system, Q\ ~ 0.32, while for a disordered state, 
Q\ = 0; however, as we will see below, = does 
not necessarily imply a disordered state. Figure 3 shows 
schematically the behaviour of {Q{k)) as a function of k 
for various values of L. For fixed wave-vector ^ 0, we 
see that {Q{k)) ^ as L ^ 00. To study the fc ^ 
limit, we fix the mode number to (e.g. to = 1 or to = 2) 
and monitor = (Qm) = {Q{2Tnn/L)). We see that 
Qh Q2 6tc. approach a finite limit as L — > 00. Numer- 
ical simulations reveal that for particles sliding down a 
fluctuating EW surface, Q* ~ 0.18 and Q2 ~ 0.09, while 
Ql ~ 0.16 and Q2 ^ 0.08 for the case of driving by a 
KPZ surface. 

A crucial point about the 1ow-to Fourier modes Qm 
is that they have broad distributions in the thermody- 
namic limit. In time, the state shows strong fluctua- 



tions in the type and amount of ordering, within a sub- 
space of ordered states. This leads us to call this phe- 
nomenon fluctuation-dominated phase ordering (FDPO). 
It is quite unlike normal ordering in equilibrium phase 
transitions, where the order parameters, such as the 
density in a liquid-vapour system or the magnetization 
in a ferromagnet, have well-defined values at all times 
in large systems. The order parameter distributions 
Prob{p) and Prob{m) in these cases consist of two delta 
functions at values characteristic of each phase, i.e. at 
well-defined densities pv and pi for the fluid, or at mag- 
netizations ±TOo for the magnet. By contrast, in our 
system, Prob{Qi), Prob{Q2) etc. approach well-defined 
broad distributions in the thermodynamic limit. 



{Q(k)) 




k — 2irm/L 



FIG. 3: Plots of {Q{k)) versus k are shown for various lattice 
sizes Li < 1/2 < L3 < L4. As L increases, for each fixed low 
value of m, — (Qm) approaches a limiting value whose 
significance is discussed in the text. 



The values of Q*, ' ' ' characterize the sort of clus- 
tering that occurs in the system. For instance, Qi is 
largest in configurations with a single dense cluster of 
particles which extends across half the system. On the 
other hand, Q2 is largest in configurations with two sepa- 
rated dense clusters. The time series for Qi shows that it 
makes large excursions around its mean value, occasion- 
ally reaching high values ~ 0.32, and occasionally low val- 
ues ~ [l3| • It is observed that a dip in the value of Qi is 
usually accompanied by a rise in the value of Q2. On the 
rare occasions when Qi and Q2 are both small, picks 
up. These observations are consistent with strong fiuc- 
tuations of the ordered structure: a single high-density 
region is most likely, but it occasionally breaks into two, 
and more infrequently into three such regions, and so 
on — never, however, lapsing into a disordered state. 
FDPO thus involves the system circulating within an at- 
tractor of states, each with macroscopic order. In view 
of the strong anticorrelations between Qi,Q2,Q3,' ■ the 
best characterization of the state would be through the 
joint probability distribution Prob(Qi, Q2, ■ ■ ■) but this 
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remains to be studied systematically. 

Evidence for this picture comes from monitoring the 
length imax (t) of the longest cluster in a system of size 
L [1^. In a disordered state, (imax) scales as logL. A 
numerical study of the probability distribution of imax 
shows that {imax) increases as L'^ where the exponent 
ranges from 0.6 to 0.9 for different sorts of surface dynam- 
ics, and the full distribution scales with (imax), implying 
that the system does not reach a disordered state. 



V. SLIDING NONINTERACTING PARTICLES 

The related problem of noninteracting particles ad- 
vected by the Kuramoto-Shivashinsky equation was stud- 
ied by Bohr and Pikovsky while Chin [H] studied 
advection by a KPZ field. In both these studies, no noise 
acts on the particles, in which case clusters do not break 
up. Besides studying cluster coalescence in time, they 
monitored the rms displacement of a tagged particle, and 
found that it grows as t^^^ where z = 3/2 is the dynamic 
exponent of the 1-d KPZ model. Drossel and Kardar [2^ 
studied the rms displacement in the presence of noise, 
and found the same behaviour, while Gopalakrishnan 
studied the same quantity for driving by an EW surface. 
Further, Drossel and Kadar studied the density-density 
correlation function in the cases of particles falling along 
(advection) and against (anti-advection) the local slope. 
In the latter case they concluded [2^ that correlations 
show a power law decay r~^) with increasing separa- 
tion r. We will discuss this further below in the light of 
our results. 



they block evolution on the sites at which they reside. In 
Section 5 below we will see that this form of the corre- 
lation function is found also in the adiabatic limit of a 
quenched interface. 

In the case of KPZ anti-advection, simulations show 
that Eqs. 9 and 10 hold, with ^ and v ~ 0.31. The 
important point is that C is a function of r/L and not 
r alone. Thus the power law r^'^ has an L-dependent 
prefactor, unlike the power law behaviour at a critical 

P'nt. In this respect, our result differs from that of 
M- 

In the case of EW evolution, Eqs. 9 and 10 are found 
to hold again, with /i = and v ~ 2/3. 

The probability that a given site holds n particles 
L) also assumes the scaling form 

P(n,L)^L-2^/(-^) (11) 

where 

I{V)^V-' as 2/^0. (12) 

Numerical simulations for the case of KPZ advection 
yield 5 ~ 1, 7 ~ 1.15; the numerical results are also 
consistent with 7 = 1 with logarithmic corrections. For 
the case of KPZ anti-advection, we find 5 ~ 1/3, 7 ~ 1.7 
while for EW surface dynamics we find 5 ~ 0.68, 7 ~ 1.5. 

Similar results were obtained in two dimensions, using 
a lattice model similar to that discussed in Section 4.4. 
The two point correlation function shows a divergence 
- {r/L)-"^ with V2 ~ 1.4 for KPZ advection, 1^2 ^ 0.5 
for KPZ antiadvection and 1^2 — 0.3 for EW dynamics. 



A. Correlation functions and number distribution 

Using Model F defined in Section 3 (see Fig. 2b) 
numerical simulations were performed to monitor the 
steady state correlation function C{r,L) — {nini^r)L 
[H, [13. Both EW and KPZ surfaces were considered 
in 1-d. Since there is no limit on the particle occupancy 
Hi at any site, the possibility arises of a much larger de- 
gree of clustering. This manifests itself in the fact that 
the correlation function is a scaling function of r and L. 
The small-argument singularity of the scaling function 
shows a divergence: 

Csir,L)=L-^^Y[^) (9) 

with 

Y{y)^y-' as 2/^0 (10) 

The divergence is strongest for the case of KPZ advec- 
tion, in which case = 3/2 and fi = 1/2. This result 
agrees with the analytical results of Derrida et al. [3l| 
for a slightly different model, which consists of two parti- 
cles which slide down slopes, but are not passive in that 



B. Adiabatic limit: the Sinai model 

Consider the adiabatic limit in which the surface is 
absolutely still, while particles evolve by performing bi- 
ased random walks, with the bias on each link set by 
the slope of the quenched surface. The problem then 
reduces to the Sinai problem of noninteracting random 
walkers moving in a random medium which itself is de- 
fined as the trajectory of a random walk. Since the sur- 
face is static, the walkers have infinite time to explore the 
landscape, and eventually reach an equilibrium state de- 
fined by a temperature T = Vhi{p2/q2)- Thus {niUi^r) 
is given by exp[— /?(/ii + ft,j+,.)]/Z({ft,fe}) for a particular 
configuration of heights {hk}- We need to further aver- 
age (niUi^r) over quenched disordered landscapes {hk}, 
a computation that was carried out by Comtet and Tex- 
ier [321 1 . Taking the scaling limit r — > 00, L — > 00 with 
r/L held fixed, we find that {niUi+r) has the scaling form 
of Eq. 2, with 

F(y) = (2n/32L)-i/2[y(i„j^)]-3/2 (^3) 
where y — r/ L. 
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Further, the probabihty P(n,L) that a site holds n 
particles can be calculated as well ^4] . It has the scaling 
form P(n,i) = A / f3'^ L"^ X (x) with x = 2n/P^L and 

Xix) ^ e-'=Ko{x)/x (14) 

where Ko{x) is a Bessel function of imaginary argument. 

Surprisingly, we find very good agreement between 
these results for an equilibrium state of particles in a dis- 
ordered environment, and those for the strongly nonequi- 
librium situation of KPZ advcction discussed in Section 
5.1. The correlation scaling function Y{y) in Eq. 13 fits 
the numerical data closely, if we set /3 ~ 4. On the other 
hand, the scaling function X{x) in Eq. 14 describes the 
probability density data for the nonequilibrium system 
provided we set P ~ 2.3. Although the driving force is 
different in the two cases - being the temperature in the 
equilibrium disordered system, and surface fluctuations 
in the nonequilibrium system - it allows particles to ex- 
plore the terrain in both cases, and reach states with 
similar though not identical characteristics. 

VI. CONCLUSION 

The principal conclusion of this paper is that in a class 
of inhomogeneously clustered states, the density-density 
correlation function is a function not of the separation 
alone, but rather the ratio of separation to system size 
(in steady state) or separation to a growing length scale 
(in the coarsening regime). This scaling behaviour is rem- 
iniscent of phase ordered states, but unlike conventional 
phase ordering, the scaling function here is singular at 
small argument. In a system of particles with mutual 
exclusion sliding down a fluctuating surface, the singu- 
larity in question is a cusp, while for noninteracting par- 
ticles which can cluster more strongly, the singularity is 
a divergence. There is a degree of universality in the 
exponents which characterize the singularities in scaling 
functions: they remain unchanged under variation of pa- 
rameters in the model, for instance, the ratio of rates of 
particle movement and surface fluctuations, though they 



do depend on the symmetry and dimension of the driving 
field. 

The other hallmark of the states under consideration 
here is the occurrence of giant fluctuations, characterized 
by the standard deviation of the order parameter grow- 
ing proportionally to the mean, as a consequence of a 
probability distribution that remains broad in the ther- 
modynamic limit. 

It is to be emphasized that the fluctuation-dominated 
states under discussion here are quite different from crit- 
ical states. In critical states, correlation functions decay 
as power-laws, and size effects appear as corrections. By 
contrast, correlation functions in our case are, even in 
leading order, functions of the separation scaled by sys- 
tem size. Moreover, fluctuations at critical points lead to 
probability distributions that narrow down in the ther- 
modynamic limit, unlike the giant fluctuations present 
here. 

The occurrence of singular scaling functions is not con- 
fined to the models studied in this paper. A study of 
nonequilibrium nematic states shows that the density- 
density correlation function shows similar scaling prop- 
erties, with a cuspy behaviour of the scaling function [15j . 
Further, a recent study of inelastically colliding particles 
for which the coefficient of restitution depends on the rel- 
ative velocity of approach of two particles, shows similar 
non-Porod scaling properties p^ . 

Finally, there is an intriguing connection between the 
properties of the strongly nonequilibrium system under 
study and an equilibrium system with quenched disorder 
which is obtained by freezing the fluctuating potential. 
We have presented evidence for the similarities of several 
properties, but a deeper understanding of the connection 
remains an open question. 
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